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We propose and analyse a GUT-inspired version of the MSSM featuring non universal sfermion 
masses boundary conditions and asymptotic 6-r Yukawa coupling unification. The particle 
spectrum of the model gives rise to extended coannihilation processes, involving the lightest 
sbottom and the tau sneutrino. We study the allowed parameter space, subject to cosmo- 
phenomenological requirements, and the details of the new coannihilations modes. Implica- 
tions for direct and indirect neutralino searches, as well as for collider physics, are also briefly 
summarized. 



1 Introduction and Motivations 



The physics at the interface of cosmology and elementary particles theory is currently undergoing 
an exciting era. On the one side, we entered an epoch of precision measurements in the field 
of astrophysics and cosmology, as demonstrated by the recent results obtained by the WMAP 
surveyl^. On the other side, the hunt for super symmetry, which would solve a number of open 
fundamental questions concerning the physics beyond the standard model, is going to be boosted 
both by accelerator - and by direct or indirect neutralino search experiments . 
As regards the theoretical and phenomenological exploration of supersymmetric models, a num- 
ber of recent experimental and observational results increasingly constrain the candidate scenar- 
ios and their parameter space. In this paper we deal with the so-called minimal supersymmetric 
extension of the standard model, in a gravity mediated SUSY breaking context!^. In these mod- 
els, the numerous a priori free low energy parameters are fixed, through RG running, by the 
values of the soft SUSY breaking masses and trilinear terms at some high energy scale Mx- 
The common practice is to assume that some underlying principle forces these high energy vari- 
ables to some universal values at Mx, thus reducing the unkowns to four numbers and one sign 
(CMSSM). 




Figure 1: The evolution of the coupling constants and of the soft SUSY breaking masses of the 5 and 10 SU{5) 
multiplets for a 5*17(5) SUSY GUT theory with Mqut smaller than the scale where scalar universality is assumed 

(here Mpianck)- 



If one further supposes that the electroweak and strong interactions ultimately unify in a single 
grand unification theory (GUT), as suggested by the running of the coupling constants, the GUT 
gauge group can dictate some relations among the Yukawa couplings of the high energy theory 
(Yukawa Unification)!^!. 

We will analyse here the phenomenology of a class of GUT-inspired models, characterized by 
the relaxation of the hypothesis of universality in the sector of the fermions scalar superpartners 
and by partial (6-t) Yukawa Unification (YU) ^. These models are motivated by a minimal 
SU{5) SUSY GUT theory in which the scale of grand unification Mqut lies below the scale 
Mx where universality is assumed^. The SU{5) RG running between Mx and Mgut induces a 
non-universality in the whole scalar sector (the gaugino sector universality is instead not spoiled 
by the runningL'J). As regards the sfermions, the multiplet structure of the GUT theory, which 
collects the particle content into the 5 (down-type quark singlet and lepton doublet) and the 
10 (up-type quark and lepton singlets and quark doublet) representations of SU{5), implies the 
following parameterization for the soft SUSY breaking masses at Mgut- 



ttlq = mfj = m\ = m^Q = ttlq (1) 
m\ = rrijj =m\= K'^tjiIq (2) 

In fig^ we sketch a typical pattern for the running of m^Q and m"^, compared above with 
the coupling constants running. Although the full GUT theory would also imply a large non- 
universality in the Higgs sector and K > 0.7, we will resort to a simplification of the model, 
and let m'j^_^ ^ = rrig in order to concentrate only on the implications of minimal non-universal 
sfermion masses (mNUSM). We also vary K between and 1 (the latter being the fully universal 
case). We therefore generically expect a lighter spectrum, with respect to the universal case, for 
the sparticles belonging to the 5 representation for K < 1. 

On top of sfermion non-universality, we inherit from the original SU{5) SUSY GUT the successful 
prediction of b-r YU. We impose exact YU at Mgut, and fix the common Yukawa coupling 
hfy = hr{MGUT) in order to get the right value for m-j-iMz), after having included running 
effects and the SUSY thresholds corrections at MsusY = ^^ii ■ "^^^ output is mb{Mz), 
which is then required to lie within the properly evolved experimental range °. 
The current upper limit on the cosmological relic density of neutralin os as cold dark matter 
candidates requires, in the CMSSM, efficient suppression mechanisms l^^^. Among these, a 
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Figure 2: The sbottom and sneutrino mass spectrum at M1/2 = 1.1 TeV, tan/3 = 38.0 and Ao = 0, for several 
values of mo, as functions of the non-universality parameter K. For the tau sneutrino, from left to right, mo = 
1350, 1650, 1950 GeV, while for the sbottom mo = 2850, 3300, 3750 GeV. 



recent uprise of interest has involved the so-called coannihilation processes , occurring when 
the next-to-lightest SUSY particle (NLSP) mass is close to the lightest, which, for cosmological 
reasons, we require to be the neutralino. The main scope of the present analysis is to shed some 
light on the possibility that the minimal proposed deviation from universality in the sfermion 
soft SUSY breaking masses could activate extended coannihilation channels, involving sparticles 
belonging to the 5 (stau, tau sneutrino, sbottom). Further, the requirement of b-T YU, and the 
fulfillment of a set of cosmo-phenomenological constraints will strongly limit the allowed tan (3 
values as well as the sparticles mass range. The implications of these bounds for accelerator 
searches and for direct and indirect SUSY CDM detection will also be briefly discussed. 



2 The Particle Spectrum 

The large values of the Yukawa couplings of the third generation fermions, owing to RG running, 
yield lighter superpartners than the ones of the other two generations. Besides the widely 
discussed case of the staul^, the role of NLSP, in presence of non-universal sfermion boundary 
conditions, can be played by both the sbottom and the tau sneutrinJ^. In order to quantify this 
statement, we plot in figElthe masses of the two sparticles at various values of ttiq as functions 
of the parameter K. The features shown in figUl are traced back to the generic form of the 
approximate solutions to the RG equations at the electroweak scale, which can be parameterized 
as 

~ moJ ao^i + bofiK"^, (3) 




Figure 3: Mass isolevel curves for the tau sneutrino and for the sbottom in the (if, mo/Mi/2) plane. The solid 
red and black lines indicate respectively ma^ = and mg = 0. The splitting between the isolevel curves is 10% 
of the neutralino mass m^ ~ 500 GeV; tan 13 = 38.0, Aq = Q, < Q. 
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Figure 4: (a) the neutralino relic density Q<^h^ and (6) the inclusive branching ratio BR(b sj) for three 
benchmark K cases, corresponding respectively to NLSP = f, b, i>r, at tan/3 = 36.0 and fixed Anlsp = 10%. 

where a and b are, to a good approximation, functions of the ratio {mQ/Mi/2) of tan (3, and 
weakly depend on the trihnear couphng ^o^- In ah cases we notice that for a given value of 
rriQ there exists a corresponding range for K where the sbottom or the tau sneutrino are quasi- 
degenerate with the neutralino " (whose mass roughly depends only on M1/2 and is therefore 
insensitive to K). The mass isolevel curves plotted in figEl highlight the plausible coannihilation 
regions produced by the particle spectrum of NUSM models (see sec|lj). The shape of figOlis 
qualitatively unchanged by the variation of any SUSY parameter, slowly moving to the right 
when increasing tan/?!^. We then conclude that we expect extended sfermion coannihilations 
for K < 0.5 and for relatively large values of niQ > 1.5 Mi/2- 

3 Cosmo-Phenomenological Constraints 

Supersymmetric models with conserved R parity are commonly required to fulfill two categories 
of phenomenological constraints. On the one side, since these models produce natural candi- 
dates for the CDM content of the Universe, in our case the neutralino, the corresponding relic 
abundance as well as the direct or indirect cosmic detection rates must be compatible with the 
current observational and experimental data. On the other side, many precision accelerator 
measurements can be used to test the viability of a given supersymmetric scenario. Finally, 
the extra theoretical constraint of b-T YU forces the mass of the bottom quark to values to be 
compared with the current experimental and lattice data. Details on the numerical calculations 
and on the statistical analysis we adopted in this work can be found in ref.^^. 
Neutralino relic density and detection. We impose to the model the current determina- 
tion of the upper limit on the CDM content of the Universe, coming from the analysis of the 
WMAP results^ namely Ocdm/i^ = 0.1126+o;oo9} . The lower limit can be evaded under the 
hypothesis that other components contribute to the CDM besides the neutralino. The numer ical 
computation of the neutralino relic density is performed with the updated micrOMEGAs codeEIl, 
which includes thermally averaged exact tree-level cross-sections of all possible (co-) annihilation 
processes, an appropriate treatment of poles and the one-loop QCD corrections to the Higgs 
coupling with the fermions. Fig^] (a) shows our results for three benchmark cases, respectively 



^As regards the question of the fine-tuning of the parameter K, see the discussion in sec. 6. 4 of ref. 




at K = 0.2, 0.4 and 0.8. In the latter case the coannihilation region overlaps the A-pole rapid 
direct annihilation funnel, extending the allowed niy^ range up to ~ 700 GeV. The other two 
cases, instead, show the typical behavior oc m|, though suppressed by coannihilation phe- 
nomena. 

As for the direct and indirect neutralino searches, we verified that in any case the detection 
rate estimates for mNUSM models lie below the current experimental sensitivity*^. We find 
nonetheless that t-channel exchanges of lighter (for K <^ 1) sfermions generate higher detection 
rates for small values of K, falling within reach of future WIMP search experiments^. 
Accelerator constraints. The most stringent bounds coming from accelerator data are found 
to be, in the present work, the inclusive BR (b — > 57) (see figEl frame (6)) and the lower 
limit on the lightest CP-even Higgs mass^^. They both give lower bounds on the mass of the 
neutralino m^. As regards the SUSY contributions to the muon anomalous magnetic moment 
6a^, we chose to resort to the so-called super- conservative approach^^, and to take as a constraint 
the extended 5-cr bound derived from the latest experimental and theoretical results^^on a^. 
The b quark mass. In the context of b-T YU, the SUSY corrections to the mass of the b quark 
constrain the allowed tan /3 range as well as the sign oi fi. In fact, the sign of these corrections is 
given by the sign of /i, and is proportional to tan (3; since the tree level value of nib, fixed by the 
common b-T Yukawa coupling determined by nir, is typically higher than the upper experimental 
bound on nib, the corrections are required to be negative (hence < 0) and large (hence large 
tan P) . In mNUSM mo delsl^, the adopted top-down YU excludes /x > and forces, in the /i < 
case, 31 < tan (3 < 45. 



4 Extended Sfermion Coannihilations 



Minimal NUSM produce a new coannihilation branch at K < 0.5 (see figEl, where the coan- 
nihilating NLSP is the tau sneutrino at low 1.5 < (mo/Mi/2) ^ 3, and the lightest sbottom, 
for (mo/Afi/2) ^ 3. We study in fig O the new extended sfermion coannihilations at j4o = 
and at an intermediate value of tan/3 = 38. Frame (a) refers to the tau sneutrino NLSP, at 
K = 0.1, while frame (6) to the sbottom, at K = 0.35. We parameterize the remaining vari- 
ables of the model, mo and M1/2, through niy and Anlsp = — ^^^^ via RG running, 

in fact, (mo,Mi/2) uniquely determines (m^,A7VL5p), and vice-versa. The NLSP mass range 
where ^l^ < 0.13, at a fixed given m^, is as large as allowed by the efficiency of coannihilation 



processes. We clearly notice that sbottom-neutralino coannihilations (frame (6)), which involve 
strong interaction processes, suppress the relic density much more effectively than in the x-^^t-t 
case (frame (a)). We also find that x-b and b-b (co-) annihilations evolve, through strong inter- 
actions, into few final SM states (respectively g-b and g-g or 6-6*) and largely dominate over x-X 
annihilations, while triple x-i^t-t coannihilations present a wider and more complex pattern of 
final states (see ref. In the left part of the figures we show the bounds dictated by the 95% 
C.L. exclusion regions of various accelerator constraints. Scanning the values of tan/3 we find 
that, for tan/3 > 45 and tan/3 < 31, the cosmologically allowed region squeezes until no points 
simultaneously fulfill 6-t YU and the accelerator constraints. 



5 Summary and Conclusions 

We showed that in a SU (5) SUSY GUT with = Km\o scalar soft SUSY masses new extended 
sfermion coannihilation modes take place. These involve unusual NLSP coannihilating particles, 
such as the tau sneutrino and the sbottom. The models, to which we apply the constraint 
of b-T YU, as inherited from the inspiring SU{5) theory, fulfill the present phenomenological 
accelerator bounds as well as the recent cosmological relic density required for the neutralino as 
cold dark matter constituent. 

As regards the detectability of the models under scrutiny, we findl^ that future direct WIMP 
search experiment^ will be sensitive to part of the highlighted coannihilation regions, namely at 
low m^. Moreover, as far as the CERN LHC pp collider reach is concerned^, we find that most of 
the tau-sneutrino coannihilation region will be within reach, while the sbottom coannihilation 
region (where tuq > 3M]^/2 ^^^id a heavier SUSY spectrum is cosmologically allowed), could 
escape detectability if m^^ > 500 GeV. 
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